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U Abstract 

Q^ We report the results of a recent study of young star clusters (YSCs) in luminous blue 

JL compact galaxies (BCGs). The age distributions of the YSCs suggest that the starburst 

2 episode in Haro 11, ESQ 185-IG13, and Mrk 930 started not more than 30-40 Myr ago. 

j/j A peak of cluster formation only 3-4 Myr old is observed, unveiling a unique sample of 

C^ clusters still partially embedded. A considerable fraction of clusters (30 - 50 %), mainly 

younger than 10 Myr, shows an observed flux excess between 0.8 and 2.2 /im. This so-called 

,-H near-infrared (NIR) excess is impossible to reproduce even with the most recent spectral 

I> synthesis models (that include a self-consistent treatment of the photoionized gas) . We have 

^^ used these YSCs to probe the very earlv evolution phase of star clusters. In all the three 

^pT host galaxies, the analysis is limited to the optically brightest objects, i.e., systems that 

^— V are only partially embedded by their natal cocoons (since deeply embedded clusters are 

probably too faint to be detected). We discuss possible explanations for this NIR excess, 

in the context of IR studies of both cxtragalactic young star clusters and resolved massive 

star forming regions in the Milky Way and in the nearby Magellanic Clouds. 



1 Introduction 

Young star clusters (YSCs) are considered a powerful tool to study the star formation history 
of their hosts. However, many assumptions and constraints on the evolutionary properties 
are needed in order to compare YSCs properties and the bulk of the host stellar population. 
In nearby galaxies, it is possible to study both the resolved stellar population and cluster 



population (e.g. Silva- Villa & Larsen 2011). In farther galaxies, however, the study of YSCs 
is conducted with an analysis of their integrated light at different wavelengths under the 
assumption of an instantaneous burst, i.e. single stellar population. 

The luminous (M^ < —17.0, corresponding to total stellar masses of ~ ]^q9-io Mq) 
blue compact galaxies (BCGs) show clear signatures of interactions and/o r mergers, and the 



numerous observed massive YSCs are likely the result of these encounters (Ostlin et al. 2003 
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Adamo et al. 2010a 2011a|b ). The very bright ultraviolet and optical luminosities of these 
systems suggest rather low dust and metallicity content. Spectra dominated by emission lines 
clearly demonstrate that BCGs are undergoing a starburst episode. The youth of the burst 
episode is also observed in the recovered age distribution of the star clusters, which shows a 
peak of cluster formation younger than 5 Myr. 

The analysis of the young cluster populations in BCGs is quite challenging due to the 
rapid evolution a cluster experiences during the first 10 Myr (still partly in an embedded 
phase). Moreover, this analysis is based on the integrated luminosities of the clusters, which 
appear unresolved at the distance of the targets. Observations of resolved newly born star 
clusters in the Milky Way and in the Magellanic Clouds reveal that these are quite complex 
systems. A cluster forms in a hierarchical medium, from the fragmentation and collapse of 



giant molecular cloud complexes (Elmegreen 2010). This implies that stars form in a frac 



tal distribution and only in some cloud cores are their number densities is high enough to 



eventually form a bound cluster (Bressert et al. 



namical studies. IGieles & Portegies Zwart (2011 ) observed a continuous distribution between 



2010). This picture is confirmed also by dy- 



associations (unbound) and clusters (bound) at very young ages, and a clear distinction at 
older stages, when the crossing time of unbound systems exceeds their stellar age. 

In a newly born cluster, the massive and short-lived stars, rapidly reach the main 
sequence and produce strong winds and UV radiation, which ionizes the intracluster gas and 
create bubbles and shells. These Hii regions surround the optically bright core of stars and 
significantly contribute to the integrated fluxes. However, a large fraction of the stars is still 
accreting material from their dusty disks (young stellar objects, YSOs) or contracting (in the 
pre-main sequence phase, PMS). Moreover, the edges of the clusters are places for triggered 
( Elmegreen|2011 ) and progressive star formation (e.g. Carlson et al.|2007 ), which can explain 
the observed age spread in the star forming regions (e.g. Bik et al. 2010). 



2 Evidence for the NIR excess 



A significant fraction of young star clusters in Haro 11 (Adamo et al. 2010a), ESO 185-IG13 



(Adamo et al. 2011a), and Mrk930 (Adamo et al. 2011b) shows a clear signature of a flux 



excess at near-IR wavebands. The models used (Zackrisson et al. 2001 Adamo et al. 2010b) 



include a self-consistent treatment of the the photoionized gas, important during the first 
few Myrs of the cluster evolution. However, these models are not able to reproduce the NIR 
observed fluxes of the clusters. In Fig. [T} we show two representative cases of cluster spectral 
energy distributions (SEDs). A large fraction of analysed clusters in the 3 galaxies have 
regular SEDs, easily fitted by our models (see left panel,Fig. [I]). However, a considerable 
number of clusters have SEDs similar to cluster # 43 in the right panel of Fig. [T] Three 
different sets of fits have been performed for each cluster (see Adamo et al.|[2010a , |2011a for 
details) and depending on the outcomes, clusters have been divided into 3 samples: 1) regular 
SEDs have been fitted from UV to IR; 2) if the source presented an excess at A > 1.0 //m, 
the fit was performed excluding the IR data (so called IR excess); 3) if the excess affected 
also / band (~ 8000 A), then only the data from UV to optical (A < 8000 A) were included 
in the fit (referred to as NIR excess). The exclusion of the data points with an excess was 
necessary in order to reproduce the observed UV and optical SEDs of the clusters. The two 
observed fluxes at 0.5 and 0.6 nm of cluster # 43 are produced, respectively, by a narrow 
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Figure 1: An example of SED analysis of 2 clusters (the name of the host in indicated in the 
inset). The two objects are representative of different cases: on the left a young cluster with a 
good fit, i.e. regular SED; on the right, a cluster affected by an excess in flux at wavelengths 
longer than 8000 A. The filled black points indicated the observed fluxes for each cluster. 
The integrated model fluxes are labelled with open squares. They sit above the best fitting 
spectrum overplotted with a solid line. Red vertical lines indicate which data points have 
been excluded from the fit. The age and mass of the cluster are shown. 



filter (F550M) centred on the continuum and a filter (F606W) which transmits Ha. A jump 
between these two filters can be used as an age indicator (only clusters younger than 10 Myr 
have Ha in emission). Clearly, # 43 is a very young cluster. However, if the NIR data points 
(including / band) are included in the fit, the age obtained for this cluster is 45 Myr (see 
Adamo et al. 2010a for details). In Table III we give the fractions of clusters with a regular 
SED or NIR excess. Color-color diagrams of the whole Mrk930 cluster population (see 



Table 1: Fraction of clusters in each SED sample, of the three galaxies. 



targets 


re 


igular SEDs 

% 


IR 


excess 

% 


NIR excess 

% 


Haro 11 

ESO 185 
Mrk930 




44 
68 
62 




14 
10 
12 


42 
22 
26 




te 


analysis) are 


shown in 




Adamo et al. 201 lb| for s comple 


Fig. 2 The R 



I color is used 

as reference. Generally, a negative R — I < color indicates ages younger than 10 Myr and 
is produced by a strong nebular contribution to the R filter, which includes the bright Ha 
line. However, this assumption is not always valid due to the flux excess in the / band. The 
clusters, in the diagrams, have different symbols depending of their observed SEDs: black 
dots indicate normal SEDs, red triangles clusters with an excess at A > 1.0 /im (hereafter, 
IR excess), and blue diamonds SEDs which deviate at A > 0.8 ^m (hereafter, NIR excess). 
We notice that clusters with a NIR excess have a, R — I color between 0.2 and 1.2 mag redder 
than the prediction made by the best-fitting SED model. The R — I color of the "blue" 
cluster-diamonds is such that, if overlooked, causes age (and mass) overestimates, affecting 
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the results of the optical based cluster analysis. The UV — R (left panel) color shows that 
the clusters detected in the FUV are young, at least younger than 30-40 Myr. Many of the 
clusters with an NIR excess are located in an area where R — I > 0, e.g., corresponding to 
ages older than 10 Myr. However, their UV color is compatible with being a few Myrs old. 
The FUV band is sensitive to the reddening. Looking at the color-color diagram, one can see 
that clusters detected in the FUV have extinctions Ay < 1.0 ( e.g., the arrow in the plot). 
Therefore, we consider these very young FUV bright clusters as systems which have already 
gone through the deeply embedded phase. The inclusion of the IR color (right panel), clearly 
shows that clusters with a flux excess (clearly young) at the redder wavelengths are mainly 
located in an area of the diagram where the R — I > and H — R > 1.0 mag, e.g. apparently 
older than 1 Gyr. The IR color of the clusters with a red excess suggests that the extinction 
in these objects should be much higher than the one predicted by the optical and UV colors. 
In other words, at the NIR wavelengths it is possible that we are probing a different, more 
deeply embedded stellar populations, indicating that these are very young clusters. 



3 Possible origin of the NIR excess 

3.1 The / band excess 

The I band excess has been found only in very young clusters (usually < 6 Myr). A viable 



explanation for this feature is the extended red emission (ERE, see for a review Witt & Vijh 



2004). The ERE is observed as a soft rising continuum between 0.6 — 0.9 /im. It is observed 
around galactic and extragalactic Hii regions and caused by a photoluminescence reaction 
on dust grains heated by hard UV radiation. Such energetic photons are mainly produced in 
short-lived massive stars. This could explain why the I band excess in our clusters is over 
after 6 Myr. 

3.2 The IR excess 

Several mechanisms can concur to make the flux at A > 1.0 /um higher than predicted by 
models. The distance of the galaxy and the resolution achieved - the best with the current 
accessible facilities - limit our studies to the integrated properties of these YSCs. However, 
observations of close-by resolved clusters and numerical predictions of stellar populations in 
clusters can give us an hint of the mechanisms which are likely producing the observed excess. 

Among the youngest and massive resolved star clusters, 30 Doradus (hereafter 30 Dor) 
in the Large Magellanic Cloud (LMC), represents the best reference to understand what a 
recently born star cluster looks like. 30 Dor is the central region of the extended Tarantula 
nebula. Multiwavelength studies of this regions have dissected the different components of 



the complex 30 Dor environment. Walborn &: Blades (1997) identified five different stellar 



populations in the region: the bright core early 0-type stars which are part of the compact 
star cluster R136; in the north and west region embedded massive YSOs; 3 more evolved 
stellar population groups in the southern and 1.0' away in the western region. The R136 
cluster has a mass of ~ 10^ Mq and is 3 Myr old. This nuclear region (< 3 pc) is dust and 
gas-free. However, the aperture radius we are using in our analysis is much larger (radius of 
~ 36 pc) than the size of R136. Our apertures are comparable to the size of the image of 
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Figure 2: Color-color diagrams of the cluster population in Mrk930 (Adamo et al. 2011b). 
Different filter combinations are compared to the R — I color (F606W-F814W). Left: UV — R 
(F140LP-F606W); Right: R- H (F606W-F160W). In each panel, the ZOl evolutionary 
tracks are plotted as a solid black line. Where predictions for the used filters were available 
we included the M08 tracks (dashed lines) as well. Ages are labelled along the tracks. The 
black dots are clusters with regular SEDs (from UV to IR) . The red triangles are cluster with 
an excess at A > 1.0 /um. The blue diamonds shows objects with an excess starting longword 
A > 0.8 /im. The black arrows indicate an visual extinction of Ay = 1. Mean errors are 
included. A solid red line divides the plots in two regions: R — I > (older than 10 Myr) 
and R — I < (younger than 10 Myr). 



30 Dor showed in Figure 1 of Campbell et al. (2010). This suggests that while the optical 



range is dominated by the stellar and gas emission, the IR transmits also flux from the diffuse 
dust, heated by the hard UV radiation field, the embedded YSOs formed in triggered star 
formation events at the edge of the nucleus, where most of the dense gas and dust filaments 
are located, and the low mass stars still in the PMS phase. 

In the literature, studied cases of IR excess in young embedded or partially embedded 
unresolved extragalactic clusters have explained the red excess as due to an important con- 
tribution by YSOs ( Fernandez-Qntiveros et al.||2009 ), or hot dust ( Cresci et al.|2005 ). Likely, 
the same mechanism is causing the excess in young star clusters in BCGs. 

After several Myr this complex phase is over, so it cannot explain why we still observe 
objects with an IR excess at older ages. For these evolved clusters a possible source of excess 
can be an important contribution from red super giants (RSGs). Models usually predict the 
number of RSGs, assuming that the stars in a cluster fully populate the underlying initial 
mass function (IMF). However, this assumption is not valid, if the cluster is less massive than 
10^ Mq, and causes important variations for cluster masses below 5 x 10^ M© (Silva- Villa & 



Larsen|2011 ). Moreover, it has been observed that in metal-poor environments the number of 
observed RSGs tends to be higher than the predicted one from the ratio of blue versus RSGs. 
Therefore, both effects would be observed as a rise in the NIR integrated flux of an unresolved 
cluster, which our current models cannot fully account for. It is not clear, however, why we 
do not see any mass dependence between the excess in H band and the mass of the cluster, 
which could support the stochasticity scenario, or why our models cannot predict the correct 
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number of RSGs only for some of the clusters. NIR spectroscopy is needed to test these 
scenarios. 

Another possible explanation is that a second stellar population is forming in the clus- 
ters. If a second population is forming now, then we expect to detect it mainly in the NIR. 
Multiple stellar populations have been detected in globular clusters and, recently, even in 



young star clusters (Larsen et al. 2011). Are we observing the formation of a new stellar 
population in these massive YSCs? 

Finally, a bottom heavy IMF, i.e. a much steeper slope at low mass ranges (higher 
number of low mass stars), could also explain the excess. However, studies of IMF variations 



in massive star clusters (see Bastian et al. 2010) have not found any confirmed case of such 



type of cluster. Moreover, it is not clear how, in the same galaxy, a fraction of clusters could 
form with a different IMF. 
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